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ABSTRACT 

We present new empirical estimates of the AV^^^ parameter for 15 Galactic 
globular clusters (GGCs) using accurate and homogeneous ground-based opti- 
cal data. Together with similar evaluations available in the literature, we ended 
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up with a sample of 62 GGCs covering a very broad range in metal content (- 
2.16<[M/H]<-0.58 dex). Adopting the homogeneous metallicity scale provided 
either by Kraft & Ivans (2004) or by Carretta et al. (2009), we found that the 
observed AVj^^^ parameters are larger than predicted. In the mctal-poor regime 
([M/H]<-1.7, -1.6 dex) 40% of GCs show discrepancies of 2a (^0.40 mag) or 
more. Evolutionary models that account either for a- and CNO-enhancement 
or for helium enhancement do not alleviate the discrepancy between theory and 
observations. The outcome is the same if we use the new Solar heavy-element 
mixture. The comparison between a- and CNO-enhanced evolutionary mod- 
els and observations in the Carretta et al. metallicity scale also indicates that 
observed AV^^^ parameters, in the metal-rich regime ([M/H]>0), might be sys- 
tematically smaller than predicted. 

Subject headings: globular clusters: general — stars: evolution — stars: horizontal- 
branch — stars: Population II 

1. Introduction 

Globular clusters (GCs) possess several evolutionary features that show up in the color- 
magnitude diagram (CMD) and/or in the luminosity function (LF). Along the red giant 
branch (RGB) there is a well defined bump in the differential LF and, cquivalcntly, a change 
in the slope of the cumulative LF. This occurs when the H-burning shell crosses the chem- 
ical discontinuity left behind by the deepest penetration of the convective envelope (first 
dredge-up) at the base of the RGB (Thomas 1967; Iben 1968; Renzini & Fusi Pecci 1988; 
Castellani et al. 1989). The efficiency of the H-burning shell is affected by the sharp increase 
in the hydrogen abundance and the stellar luminosity experiences a temporary drop. Stars 
thus cross the same luminosity range three times, and a bump/overdensity appears in the 
differential LF of RGB stars. 

The RGB bump is interesting because: i)- the luminosity and shape provide robust 
constraints on the chemical profile inside RG structures and demark the maximum downward 
penetration of the convective envelope during the first dredge-up (Bcrgbusch & VandcnBcrg 
2001; Cassisi et al. 2002); ii)- the luminosity peak can also be used to estimate cither 
the cluster metallicity (Desidera et al. 1998) or the cluster distance. However, the bump 
luminosity is affected by both theoretical (treatment of the convective transport in stellar 
envelopes) and observational (distance, reddening) uncertainties. To overcome the latter, 
Fusi Pecci et al. 1990, hereinafter FP90) suggested using the AVH^^^V{fMmp) - V{HB) 
parameter, i.e., the difference in apparent visual magnitude between the RGB bump and 
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the Horizontal Branch (HB) at the luminosity level of RR-Lyrae instability strip. Using 
new estimates of this parameter for 11 GCs, FP90 found that the observed AVj!^^^ values 
were 0.4 mag fainter than predicted. This gave rise to a series of theoretical (Alongi et al. 
1991; Bono & Castellani 1992; Bergbusch & VandenBerg 2001) and observational (Alves & 
Sarajedini 1999; Ferraro et al. 1999) investigations. 

Cassisi & Salaris (1997) gave a new spin to this problem, finding good agreement between 
theory and observations once the global metallicity of individual clusters was taken into 
account. This was confirmed by Zoccali et al. (1999) using a large sample (28) of GCs 
covering a wide metallicity range (— 2.1<[Fe/H]<— 0.3 dex). The comparison between the 
evolutionary lifetimes during the crossing of the H-discontinuity and star counts across the 
RGB bump led also to good agreement between theory and observations (Bono et al. 2001). 
More recently, Riello et al. (2003, hereinafter R03) found that predicted AV^^^ values agree 
quite well with observations for cluster ages of 12±4 Gyr. This investigation used a sample 
of 54 GCs, but only three metal-poor ([M/H]<-1.75 dex) GCs were included. Note that 
identification of the RGB bump in metal-poor GCs is more difficult since it is brighter than 
the HB, where evolution along the RGB becomes faster and the stellar sample, consequently, 
becomes smaller. 

Here we augment the previous observations by investigating AV^^^ in a large sample of 
GCs with metallicities ranging from —2.43 to —0.70 dex. In particular, we focus our attention 
on different sources of possible systematic errors, namely metallicity scale, a elements, CNO 
abundances and helium content. 



2. Empirical and theoretical frameworks 

From the database maintained by Stetson (2000) we selected 15 GCs with low reddening 
{E{B-V) < 0.10, Harris 1996) and a wide range in metalhcity (-2.43<[Fe/H]<-0.70, see 
data listed in Table 1). The B and V data adopted in this investigation come from original 
and archival observations which have been collected, reduced, and calibrated by PBS in an 
ongoing effort to provide homogeneous photometry on the Landolt (1992) photometric sys- 
tem for a significant fraction of GCs. More details concerning individual fields can be found at 



the following URL: http://www4.cadc-ccda.hia-iha.nrc-cnrc.gc.ca/community/STETSON/standards/ 



or by direct communication with PBS. 

To find the RGB bump we represented each RG star detected in the B, y-bands by a 
Gaussian kernel with a a equal to the photometric uncertainty. The differential RGB LF was 
formed by summing the individual Gaussians and the position of the bump was estimated by 



-4- 



fitting the resulting peak with a Gaussian. Data plotted in Fig. 2 show that the photometric 
error in color, added in quadrature, is on average of the order of 0.01 mag from bright RGs 
to stars fainter than the main sequence turn-off. Therefore, the uncertainty in the measured 
position was assumed, as a generous estimate, equal to the sigma of the fitting Gaussian. 
The fit was performed using an interactive program that performs analytical fits to the main 
peaks of a given distribution (Calamida et al. 2009). On the basis of the residuals between 
the RGB LF and the analytical fit, the software allows us to remove/insert new Gaussian 
components manually. 
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Fig. 1. — Top - RGB luminosity function for a metal-poor (M92, [Fe/H]=-2.38 dex, left) 
and a metal-intermediate (NGC 5904, [Fe/H]=-1.32 dex, right) GC. The star counts (A^^) 
are normalized to the number of stars in the RGB bump. The arrows mark the position of 
the RGB bump. Bottom - (V, B-V) CMDs of the two GCs plotted in the top panels. The 
thin solid lines show the HB level, while the arrows indicate the RGB bump. 



The HB level for GCs with either an intermediate or a red HB morphology was estimated 
using the LF of HB stars in the fiat region. Specifically, the ZAHB level was estimated at 
3(7 fainter than the peak of the LF (see bottom right panel of Fig. 1). Note that the 3a 
values range from ~0.09 in metal-poor clusters to ~0.18 mag in metal-rich clusters. This 
approach may be compared to the procedure adopted by Zoccali et al. (1999) and by Riello 
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Fig. 2. — Intrinsic photometric error in color (B-V) as a function of V magnitude for M92 
(top) and NGC5904 (bottom). In each cluster the outliers at ^ ~ 15 are RR Lyrae variable 
stars. 

et al. (2003) for GCs with a red HB morphology (metal-rich). They fit the lower envelope 
of the HB stellar distribution and the ZAHB level was fixed at 3a above the lower envelope. 
The former approach is less affected by accidental measuring errors, since it relies on the 
peak of the LF. On the other hand, for GGCs with a blue HB morphology (typically metal- 
poor and metal-intermediate GCs with only a few red HB stars) we used a template cluster 
with the same metaUicity and a weU populated HB. The HB level was estimated by shifting 
the template cluster in magnitude and in color (Buonanno et al. 1986; Ferraro et al. 1992; 
Zoccah et al. 1999). The observational errors for the CMDs of Fig. 1 are presented in Fig. 2. 
The uncertainty of AV^^^ was estimated by adding in quadrature the standard errors of 
the bump and HB luminosities (see Table 1). To increase the sample size we also included 
the AV^'^*' of R03. However, their AV^'^*' parameters were based on the F555W-band of 
WFPC2@HST (see their Table 1). We have nine clusters in common with R03 and we found 
a mean difference of 0.06 ± 0.03, i.e., a factor of two smaller than the typical uncertainty. 
We applied this shift in magnitude to the R03 estimates; our results are independent of this 
magnitude correction. For the clusters in common with R03, we adopted our AV^^^ values, 
since the number of RGB stars across the bump region is larger. We ended up with a sample 
of 62 GCs. 

To compare theory and observations we adopted different metaUicity scales. 1) The 
metaUicity scale by Carretta et al. (2009, hereinafter C09), based on iron measurements of 
Fel Unes from intermediate resolution {R Ki 20,000) spectra collected with GIRAFFE at the 
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Very Large Telescope (VLT) and from high-resolution spectra collected with UVES at VLT 
[R f» 40,000). They observed 19 calibrating GGCs and for each cluster they collected ~100 
GIRAFFE and ~10 UVES spectra. The estimated mean uncertainty of these abundances 
ranges from 0.03 (NGC 1904) to 0.08 (NGC 2808) dex. 2) The metallicity scale by Kraft 
& Ivans (2003,2004, hereinafter KI03 and KI04) based on high- resolution spectra {R ~ 
30,000, Shetrone & Keane 2000) and on MARCS atmosphere models. This scale has two key 
advantages: i) iron abundances based on measurements of Fell lines, which are minimally 
affected by non-LTE effects (Thevenin & Idiart 1999); and ii) robust determinations of 
surface gravity and effective temperature. The accuracy of this scale is on average better 
than 0.1 dex. 3) Metallicities by Rutledge et al. (1997, hereinafter RHS97), derived from 
homogeneous and accurate measurements of the calcium triplet (CaT), referred to both the 
Zinn & West (1984, hereinafter ZW84) and the Carretta & Gratton (1997, hereinafter CG97) 
metallicity scales. The typical precision of these metallicities is 0.10-0.15 dex (KI04). Note 
that the cluster metallicities provided by ZW84 were based on various diagnostics (G band. 
Call K, and Mgl lines), and their typical precision is ~0.15 dex. 

To estimate the global metallicity ([M/H]) of individual GCs, we adopted the relation 
by Salaris et al. (1993) and assumed [a/Fe]=-|-0.4 dex for all GCs. This assumption is based 
on spectroscopic measurements of a-elements in GCs (see Fig. 4 in Gratton et al. 2004). 
The idea of a constant [a/Fe] as a function of the metal content relies only on three metal- 
rich GGCs (NGC 6528, NGC 6553, Liller 1), while field stars show a steady decrease in 
a-element abundances in the approach to solar chemical composition. However, we adopted 
a constant [a/Fe] ratio, because field and cluster stars show different abundance patterns 
(see §3). Among our 62 GCs, metallicity estimates are available for 37 objects in RIIS97 
and 60 in KI04, while the entire sample is included in the C09 metallicity compilation. 

To compare theory with observations we adopted the a-enhanced isochrones and Zero 
Age Horizontal Branches (ZAHBs) available in the BaSTI databas^. We adopted isochrones 
that neglect atomic diffusion (Pietrinferni et al. 2006). The global metallicity ranges from 
[M/H]=-2.27 to +0.06 dex, the primordial helium content by mass is 1^=0.245, and the 
helium-to-metal enrichment ratio, AY/AZ, is 1.4, while the stellar mass of the progeni- 
tors (Main Sequence Turn Off, MSTO) ranges from 0.78 ([M/H]=-2.27) to 0.90 M© 
([M/H] =+0.06). Predicted AV^^^ values were estimated assuming a GC age of 12 Gyr 
(VandenBerg et al. 2006). 

To estimate the typical uncertainty in the predicted AVj!^^^ values we assumed an 
uncertainty of 0.1 dex in the heavy element abundances (Fe, a-elements), 0.02 in the helium 



Evolutionary models can be download from 'http: / /www. oa-teramo. inaf.it /BASTI 
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Fig. 3. — AV^*^^ as a function of the global metallicity [M/H]. From top to bottom the 
same data are plotted using different metallicity scales: Rutledge et al. (1997, RHS97) in 
the Zinn & West (1984, ZW84) scale; RHS97 in the Carretta & Gratton (1997, CG97) scale; 
Kraft & Ivans (2004, KI04); C09. The number of GCs in each panel is also labeled (Nc). The 
circles display AV^^^ parameters from Riello et al. (2003), while the squares show current 
estimates. The red line shows predictions (BaSTI) at fixed age (12 Gyr) for a-enhanced 
chemical compositions. The error bar plotted in the right corner of the top panel shows 
theoretical uncertainties. The GCs with difference from theory larger than 2a and 3a are 
plotted in green and in blue. 

content by mass (Y), and 1 Gyr in the cluster age. Moreover, we also included uncertainties 
in the input physics used in constructing evolutionary models (Cassisi et al. 1998; Bergbusch 
& VandenBerg 2001). 
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3. Comparison with a- and CNO-enhanced models 

Fig. 3 shows the comparison between empirical and predicted AV^j^^ values in the four 
different metallicity scales. The figure shows that the observed AV^^^ values of metal- 
poor and metal-intermediate GCs are systematically larger than predicted. This discrep- 
ancy does not depend on the adopted metallicity scale and becomes strongest in metal-poor 
([M/H]<-1.5) GCs. To be more quantitative we calculated for each GC the quadrature sum 
of theoretical and empirical uncertainties (see Table 1 and Fig. 3) . We found that the GCs 
with a discrepancy larger than 2a (green symbols in Fig. 3) are two, four and five for the 
RHS97[ZW84], the KI04 and the C09 metallicity scale. It is worth mentioning that the dif- 
ference between theory and observations is ;^0.40 mag. We also found that one GC in KI04 
shows a discrepancy larger than Sa {^0.60 mag). This suggests that metallicities based ei- 
ther on KI04 or on C09 scale show a very similar behavior when compared with observations. 
However, the bottom panel shows that several metal-rich ([M/HjRi-0.5 dex) GCs in the C09 
metallicity scale are located below the predicted loci. Note that possible evolutionary effects 
in the HB phase would imply an increase (more positive) in the value of the AV^^^ param- 
eter. Fig. 3 indicates that — if we trust the KI04 and C09 metaUicity scales — approximately 
40% of the metal-poor ([M/H]<-1.7, -1.6 dex) GCs show a discrepancy with theory that is 
at least at a 2a level. Such a discrepancy can hardly be explained as a spread in cluster age: 
the age derivative of the AV^^^^ parameter is ^ 0.03 mag per Gyr over the entire metallicity 
range (see, e.g., Cassisi & Salaris 1997; Zoccali et al. 1999), and an age spread of order 10 Gyr 
seems out of the question. Moreover, changes in mixing-length efficiency have little effect 
on the luminosity of the RGB bump: a change of ±0.1 in {ami) causes a variation of ~0.03 
mag m the AVj^^^ parameter (Cassisi & Salaris 1997). The comparison between predicted 
and empirical cluster RGBs does not support larger variations in the a^i parameter. The 
anonymous referee raised the problem that we are using a sample of AV^^^ parameters 
that are based on two different methods to estimate the ZAHB luminosity. According to the 
referee there is a big possibility of systematic errors in the ZAHB determinations. However, 
some circumstantial evidence argues against this working hypothesis. 

i)- The photometric precision of the current data and of the HST data adopted by Riello 
et al. at the ZAHB level is better than 0.010-0.015 mag over the entire cluster sample. We 
have nine GCs in common with Riello et al. that cover the entire metallicity range and the 
difference in the estimated AV^^^ parameters is smaller than 0.1 mag. 
a)- The method we adopted to estimate the ZAHB luminosity (3(7 fainter than the peak of 
HB LF) might be prone to underestimate the real ZAHB level. However, a less conservative 
estimate would imply brighter ZAHBs, and in turn larger (more positive) AV^^^ parameter 
and therefore a larger discrepancy between theory and observations. 

Hi)- We estimated the AV^^^ parameter for a good fraction (six out of eleven) of metal- 
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poor ([M/H]<-1.7, -1.6 dex) GCs. If we do not include the GCs measured by Riello et al. 

the discrepancy between theory and observations, in this metaUicity range, becomes larger. 

iv) - The majority (nine out of twelve) of the AV^^^ parameters of metal-rich ([M/H]>-0.3 
dex) GCs have been estimated by Riello et al. 

v) - To overcome subtle uncertainties in empirical estimates and theoretical predictions we 
are focusing our attention on discrepancies at the level of several tenths of a magnitude. 



- fiV;sHB=VzAHB(ouO-VzAHB(^e^^O^O St Ol. 1999) 
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Fig. 4. — Top - Difference between current Vzahb estimates and similar estimates by 
Ferraro et al. (1999) as a function of the iron abundance (C09). The dashed hnes display 
the uncertainty (±0.1 mag) we adopted for Vzahb estimates. The error bars account for 
individual uncertainties on Vzahb (sum in quadrature) and in iron abundance. Bottom - 
Same as the top, but for the iron abundaces based on the metaUicity scale by KI04. 

Finally, we also compared current estimates of Vzahb with similar estimates provided 
by Ferraro et al. (1999) to test whether our approach may also introduce a zero-point shift or 
a systematic trend with metaUicity. The quoted authors to estimate Vzahb, in a large sample 
of GGGs with non homogeneous photometry, applied a metaUicity dependent correction to 
Vhb ranging from 0.06 to 0.16 mag. We have ten GGGs in common with Ferraro ct al. and 
data plotted in Fig. 4 show no evidence of a trend with metaUicity and the two different sets 
of Vzahb estimates agree within la. 

There are several culprits that can explain the difference between current predictions and 
those provided by Gassisi & Safaris (1997). Their primordial He-content was Y—0.2S while 
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ours is y=0.245, and the estimated cluster age has been decreased from 15 to 12 Gyr. The 
input physics (opacities, equation of state) was also changed, and they adopted a scaled-Solar 
heavy-element mix rather than an a/p/ia-enhanced one. However, the difference between the 
old and the new scaled-solar AV^™^ parameter for [M/H]=-1.79 and t=15 Gyr is small 
(-0.65 vs -0.60). The difference with the new a-enhanced AV^^^ parameter, for the same 
chemical composition and an age of 12 Gyr, is also modest (—0.65 vs —0.75). Given the lim- 
ited changes between the current and previous predictions, the discrepancy between theory 
and observations in the metal-poor tail is confirmed and reinforced by the inclusion of new 
GCs in this metallicity regime. 




Fig. 5. — Same as Fig. 3, but for two different metallicity scales: Kraft & Ivans (2003, KI04, 
top) and Carretta et al. (2009, C09, bottom). The green line shows predictions (BaSTI) at 
fixed age (12 Gyr) and for a- and CNO-enhanced chemical compositions (-0.07<[M/H]<- 
1.27 dex). The blue asterisk marks the predicted AV^^^ parameter for an age of 11 Gyr and 
a metal-poor ([M/H]=-1.77 dex) a- and CNO-enhanced chemical composition. The dashed 
line shows a linear interpolation between the BaSTI predictions and the new computed 
metal-poor value. 



We note in passing that theoretical models indicate that the RGB bump luminosity 
is marginally affected by plausible assumptions concerning the smoothing of the chemical 
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discontinuity (Bono et al. 1992). Following a different approach, Alongi et al. (1991) found 
that the inclusion of nonlocal overshoot at the base of the convectivc envelope causes, in 
low-mass RGB stars, an increase of 0.4 V mag in the location of the RGB bump. However, 
this systematic shift minimally depends on the chemical composition and applies to both 
metal-poor and metal-rich stellar structures. 

Spectroscopic measurements suggest variations in the abundance pattern in GCs on 
a star-by-star basis. Together with changes in the relative abundances of CNO elements, 
well defined anti-correlations have been found between O and Na and between Mg and Al 
(Gratton, Sneden, & Carrctta 2004, and references therein). This applies not only to evolved 
and unevolved cluster stars (Briley et al. 1996; Cohen & Melendez 2005) but also to RG 
stars characterized by different thicknesses of the convective envelope (Pietrinferni et al. 
2009). This solid result and the recent evidence that field halo stars do not show, at fixed 
iron content, similar changes in (0,Na,Mg,Al) suggests the occurrence of an intra-cluster 
pollution mechanism (Denissenkov & Weiss 2004; Ventura & D'Antona 2005; Salaris et al. 
2006). Moreover, Yong et al. (2009) found that the [C-|-N-|-0/Fe] abundance ratio changes 
by a factor of four (0.6 dex) among bright RGs in the GC NGC 1851. This is consistent with 
the scenario suggested by Cassisi et al. (2008) to explain the two subgiant branches detected 
in this cluster (Milone et al. 2008) as a difference in the mixture of C+N+0 elements. 

Therefore, to investigate the discrepancy between theory and observations of AV^'^"*', 
we also considered cluster isochrones and ZAHBs constructed assuming CNO-Na abundance 
anti-correlations (Pietrinferni et al. 2009). The mixture adopted in these evolutionary mod- 
els includes a sum in CNO abundances that is a factor of two larger than in a canonical 
CK-enhanced mixture and within the observed range of GC anti-correlations (Salaris et al. 
2006). Predictions by Pietrinferni et al. (2009) are available for stars with metallicities from 
[M/H]=-1.27 {Mmsto^0.79Mq) to [M/H]=-0.07 (Mmsto^O.SSMq) dex. To constrain the 
predictions in the metal-poor regime, we used an updated version of the FRANEC evolution- 
ary code (Cariulo et al. 2004; DeglTnnocenti et al. 2008) to compute stellar isochrones and 
ZAHBs for [M/H]=-1.77 ([Fc/H]=-2.32, [a/Fe]=0.3,r=0.246, Mmsto=0.78Mq) accounting 
for both a- and CNO-enhanccmcnt. These models include atomic diffusion, which causes a 
decrease in the inferred cluster age of ~1 Gyr (Bahcall & Pinsonneault 1995; Castellani et 
al. 1997). Therefore, we estimated the AV^™^ parameter for a cluster age of 11 Gyr and 
the blue asterisk marks its position in Fig. 5. A more detailed discussion concerning the 
input physics of these models will be addressed in a forthcoming investigation (Di Cecco et 
al. 2010). Note that current predictions for CNO-enhanced models do not cover the metal- 
licity range -1.77<[M/H]<-1.27 (see dashed line in Fig 5). However, the linear interpolation 
(see dashed line in Fig. 5) is supported by the smooth change showed by more metal-rich 
structures (Pietrinferni et al. 2009). To compare theory and observations, the global metal- 
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Fig. 6. — Synthetic RGB luminosity functions for a metal-poor ([Fe/H]=-2.32, [Q;/Fe]=0.3, 
y=0.246) GC. The total number of synthetic stars is 2,500 per numerical simulation. From 
top to bottom the cluster stellar content includes fractions of 20%, 30%, 40% and 50% of 
GNO-enhanced stars. Simulations were performed assuming either the same age (11 Gyr) for 
the canonical and the GNO-enhanced sub-populations (left panels) or a difference of 1 Gyr 
between the two sub-populations (right panels). 



licity of GCs was estimated using the Safaris et al. (1993) relation with [ci;/Fe]=0.4 and 
[GNO/Fe]=0.3. The comparison indicates that a- and GNO-enhanced predictions display a 
similar discrepancy with empirical estimates in the K104 (three GGs with a difference larger 
than 2cr, and one with a difference larger than 3o") and G09 (four GGs with a difference larger 
than 2(7, and one with a difference larger than 3(t) metallicity scale (see the green lines in 
Fig. 5). The difference is mainly due to the fact that the RGB bump is more sensitive to an 
increase in CNO abundances than the ZAHB (see Fig. 1 and 3 in Safaris et al. 2008). On the 
other hand, the bottom panel of Fig. 5 shows that, if we adopt the C09 metaUicity scale, a 
group of six GCs in the metal-rich regime ([M/H]>0.0) shows systematically smaller (more 
negative) AV^^^ values than predicted. Gurrent predictions do not cover this metallicity 
range and we linearly extrapolated the trend from more metal-poor structures. Note that 
Fig. 5 shows only the KI04 and the G09 metallicity scales because they include the largest 
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fraction of GCs in our sample. The main outcome of the above comparison is that the 
observed AV^*^^ parameters in the metal-poor regime are systematically larger (more posi- 
tive), independent of the metallicity scale, than predicted by a- and CNO-enhanced models. 
A similar discrepancy, but with observed AV^^^ parameters systematically smaller than 
predicted, might be also present in the metal-rich tail ([M/H]>0.0), if we adopt the C09 
metallicity scale. 

The anonymous referee suggested a test to constrain the fraction of CNO-enhanccd 
stars that could be identified as a secondary RGB bump. Accordingly, we computed a series 
of synthetic CMDs assuming a mctal-poor chemical composition (y=0.246, [Fc/H]=-2.32, 
[a/Fe]=0.3, Pisa evolutionary code) and a cluster age of 11 Gyr. Together with this canonical 
stellar population we also included different fractions of a CNO-enhaced stellar component. 
We assumed for the CNO-enhanced population two different ages, namely 11 and 10 Gyr (see 
the left and the right panels of Fig. 6). Moreover, we have made our best efforts to match the 
synthetic CMD to the observed CMD of M92. In particular, we randomly distributed the 
same number of stars (2,500, V > mag) along the predicted RGB to mimic the Poisson 
uncertainties properly and accounting for photometric errors (ctb, cry =< 0.015 mag). Fig. 6 
shows the predicted LFs with CNO-enhanced fractions from 20% (top) to 50% (bottom). To 
define the RGB bump we adopted the same procedure as for observed GCs. Data plotted 
on the left column (coeval populations) indicate that CNO-enhanced stars show up as a 
secondary RGB bump when its fraction is equal to or larger than 30%. Note that the two 
peaks differ in magnitude by 0.1 mag, i.e. an order of magnitude larger than the typical 
photometric error. On the other hand, if we assume an age difference of 1 Gyr between the 
canonical and the CNO-enhanced population (right column) the secondary population shows 
up as a second RGB bump only when its fraction becomes of the order of 40-50%. This is 
due to the fact that the effects of a decrease in cluster age and of an increase in the global 
metallicity almost cancel out across the RGB bump. The magnitude difference between the 
peaks is 0.08 mag. This indicates that RGB bump is not a good indicator of the presence of 
a sub-population that is both CNO-enhanced and younger relative to the main population. 



4. Comparison with He-enhanced models 

Recently, accurate and deep photometry has revealed multiple stellar populations in 
several massive GCs. Together with the highly complex system u Ccntauri (Anderson 2002; 
Bcdin et al. 2004) multiple stellar sequences have been detected in GCs covering a broad 
range of metal content: NGC 2808 (D'Antona et al. 2005; Piotto et al. 2007), M54 (Siegel et 
al. 2007), NGC 1851 (Milone et al. 2008) and M4 (Marino et al. 2008). Some of these multiple 
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Fig. 7. — Same as Fig. 3, but for two different metallicity scales: Kraft & Ivans (2004, KI04, 
left) and Carretta et al. (2009, C09, right). Top - The lines display predictions (BaSTI) at 
fixed cluster age (12 Gyr), a-enhanced chemical composition and different helium contents: 
y =0.245 (red), y=0.30 (blue) and y=0.35 (green). Bottom - Same as the top, but the 
predictions refer to cluster isochrones of 10 Gyr. 

sequences {u Gen, NGG 2808) might be explained with a He-enhanced stellar population 
(Norris 2004; D'Antona et al. 2005,2008; Piotto et al. 2007). 

To evaluate the impact that the He-content has on AVj^^^ we also considered evolution- 
ary models assuming higher helium abundances. The top panels of Fig. 7 show the compari- 
son, at fixed cluster age (12 Gyr), between observed and predicted (BaSTI) AV^^^ parame- 
ters. Together with the canonical a-enhanced models (red line) the He-enhanced models are 
plotted with a blue (y=0.30; Mmsto = 0.72Mq at [M/H]=-1.27 and Mmsto=0-78Mq at 
[M/H]=-0.35)) and a green (F=0.35; Mmsto = O.68M0 at [M/H]=-0.66 and Mmsto=0-72 
Mq at [M/H] =-0.35) line, respectively. 

The He-enhanced models predict larger AVj^^^ parameters when compared with canon- 
ical models. The difference is most evident in the metal-rich tail; it is negligible in the 
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metal-intermediate regime because an increase in helium content — at fixed cluster age — 
gives brighter ZAHBs and brighter RGB bumps. Note that the He-enhanced models do not 
exist in the metal-poor regime, since these models experience either no first dredge-up, or 
only a mild one. The RGB bump vanishes together with the chemical discontinuity in the 
envelope. 

Finally, we also tested the impact of simultaneous changes in both cluster age and he- 
lium content. The bottom panels of Fig. 7 show the comparison between canonical and 
He-enhanced models for a cluster age of 10 Gyr (F=0.245, Mmsto=0.82Mq for [M/H]=- 
2.27 and Mmsto=0-9^Mq for [M/H]=+0.06; F=0.30, Mmsto=0-75Mq for [M/H]=-1.79) 
and Mmsto=0.S1Mq for [M/H]=-0.35; Y=0.S5, Mmsto=0-Q9Mq for [M/H]=-0.96 and 
Mmsto=0-74:Mq for [M/H]=-0.35). Once again the predicted AV^T" values are system- 
atically smaller than observed, regardless of the metaUicity scale. The decrease in cluster 
age expands the range in metal content within which the He-enhanced models display the 
RGB bump. However, in the metal-intermediate regime (-1.7<[M/H]<-0.8) the difference 
between canonical and He-enhanced (y=0.30) models is negligible. This is because the 
He-enhancement causes similar changes in Vhb and Vbump that cancel out in the difference. 

The anonymous referee suggested a test to constrain the impact of a sub-population of 
He-enhanced stars on the shape and width of the RGB bump. Accordingly, we computed 
several scries of evolutionary models (using the Teramo evolutionary code) at fixed global 
mctallicity, but steadily increasing the initial He abundance from the canonical value {Y ~ 
0.25) to the disappearance of the bump (sec Table 2). Fig. 8 shows the predicted width in 
visual magnitude between the faintest and the brightest point attained by the RGB bump for 
metal-poor (dashed-dotted line), metal- intermediate (dashed line) and metal- rich (solid line) 
stellar structures constructed assuming different helium contents and fixed cluster age. This 
figure shows that the width in magnitude increases by ~ a factor of five when moving from 
metal-poor to metal-rich structures. In metal-poor structures an increase in helium of ~ 10% 
causes the disappearance of the bump, since its width becomes smaller than one hundreth 
of a magnitude. The same result occurs in metal-intermediate and metal-rich structures, 
but the required increase in helium content is of the order of 30-35%. The disappearance 
of the RGB bump is a consequence of the decrease in the greatest depth of the convective 
envelope, and in turn in the reduced amount of nuclear processed material dredged-up into 
the outermost regions (Salaris & Cassisi 2008). This also means that an increase in hehum 
content has a stronger impact on the width of the RGB bump of metal-poor structures than 
of the metal-rich ones. In order to be more quantitative on this key dependence we also 
computed, for the same structures, the increase in surface helium abundance after the first 
dredge-up (see column 6 in Table 2). We found that in a metal-poor structure with F=0.35 
the amount of extra helium, i.e., the amount of helium dredged-up into the surface layers 
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Fig. 8. — Predicted width in visual magnitude of the RGB bump for metal-poor (dashed- 
dotted line), metal-intermediate (dashed line) and metal- rich (solid line) stellar structures 
computed at fixed cluster age and global metallicity, but assuming different helium contents 
(Teramo evolutionary code). 

by the sinking of the convective envelope, decreases by a factor of two when compared with 
structures constructed assuming a canonical helium content. The same outcome applies for 
metal-intermediate and metal-rich structures, but only for F=0.40. 

5. Discussion 

New electron-conduction opacities provided by Cassisi et al. (2007) have a negligible 
effect on the luminosity of the RGB bump, but affect the size of the helium core, and in 
turn the HB luminosity. The helium core mass of the new models decreases, independent of 
the metallicity, by 0.006 Mq. The difference in the V band is ~ 0.04-0.07 mag. However, 
this goes in the opposite direction, since the AV^*^"^ parameter would attain smaller (more 
negative) values. 
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Interestingly, the drift between theory and observations becomes even more evident 
when moving from the metal-intermediate to the metal-poor regime. To further constrain 
this effect, we also investigated the impact of possible uncertainties in the radiative opacities 
(Guzik et al. 2009). In particular, we constructed two new sets of isochrones and of ZAHB 
models for [M/H]=-2.27 (Z=0.0001) and [M/H]=-1.79 (Z=0.0003) by artificially increasing 
the radiative opacity by 5%. The AV^^^^ parameters based on opacity-enhanced models 
are smaller, but the difference with canonical predictions is of the order of a few hundredths 
of a magnitude. There are similar changes in Vhb and Vbump magnitude that cancel in the 
difference. 

The current discrepancy between theory and observations is little affected by the new 
solar heavy-element mixture provided by Asplund et al. (2005, hereinafter AOS). The new 
mixture has a twofold impact on GC metallicities, since it affects both the iron content 
([Fe/H]) and the relation between the iron and the global metallicity Z, i.e., the heavy- 
element abundance adopted in evolutionary computations. To quantify the difference be- 
tween the old and new solar mixtures on the iron measurements we consider the GG M92. 
According to KI04 this GG has an iron content of [Fe/H] =—2.38 dex assuming a solar iron 
abundance in number of log ei7e=7.51 (Grevesse & Noels 1993, hereinafter GN93). On the 
other hand, the iron abundance of M92 becomes [Fe/H]=— 2.32 dex, if we adopt the new 
solar iron abundance log £^6=7.45 (A05), since [Fe/H] ^05= [Fe/H] cArga + Aloge (GN93 - 
A05)= [Fe/H]GAr93 + 0.06 dex. The new Solar mixture by AOS gives global metaUicities, 
Z, that at fixed iron content are approximately 30% lower than those based on the GN93 
Solar mixture, since {Z / X)-^^^ / {Z / X)*^^^^ Ri 0.7. Given the above differences, the global 
metallicity of M92, assuming an ct-enhancement of -1-0.4, changes from Z—1.51 -10"^ (GN93) 
to Z=l. 19-10"^ (AOS). 



6. Summary and final remarks 

We estimated the AV^^^ parameter for IS GGGs using accurate and homogeneous 
optical databases. Gombining the new empirical evaluations with similar estimates provided 
by Riello et al. (2003), we ended up with a sample of 62 GGGs covering a very broad range 
in metal content (-2.16<[M/H]<-0.58 dex) and structural parameters. The comparison 
between theory and observations produced the following findings: 

a)- The observed AV^^^ parameters, using the homogeneous metallicity scales for 
GGGs provided by Kraft & Ivans (2004) and by Garrctta ct al. (2009), are larger than 
the predicted ones. In the metal-poor regime ([M/H]<-1.7, -1.6 dex) 40% of GGs show 
discrepancies of 2a (fti0.40 mag) or more. Evolutionary models that include either a- and 
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CNO-enhancement or helium enhancement do not alleviate the discrepancy between theory 
and observations. 

b ) - The above discrepancy between theory and observations is not affected by the new 
Solar heavy-element mixture. 

c) ~- The comparison between a- and CNO-enhanced evolutionary models and observa- 
tions in the Carretta et al. metaUicity scale also indicates that observed AV^^^^ parameters, 
in the metal-rich regime ([M/H]>0), might be systematically smaller than predicted. 

d) - Evolutionary models predict, as expected, a strong dependence of the RGB bump 
on the helium content. Current predictions indicate that an increase in helium of ~ 10% 
causes the disappearance of the RGB bump in metal-poor structures. The outcome is the 
same in metal-intermediate and in metal-rich structures, but the required increase in helium 
content is of the order of 30-35%. 

In this investigation we assumed a constant CNO-enhancement over the entire metal- 
licity range. However, spectroscopic evidence indicates that the CN-strength is correlated 
with Na, Al and Mg in metal-poor, with Na and Al, in metal-intermediate and with only 
Na in metal-rich GCs (Smith & Wirth 1981; Gratton et al. 2004). This implies that the 
CNO-enhancement might depend on the metal abundance (Cottrell & Da Costa 1981). 

More detailed investigations concerning the input physics adopted in the computation 
of low-mass evolutionary models is required before we can reach a firm conclusion. This 
also applies to the metaUicity scale, an in particular to the occurrence of possible systematic 
errors when moving from very metal-poor to metal-rich GCs. 

We also thank an anonymous referee for his/her comments that helped us to improve 
the readability of the manuscript. This project was supported by Monte dei Paschi di Siena 
(RL: S. Degl'Innocenti), PRIN-MIUR2007 (P.I.: G. Piotto). MZ acknowledges support 
by Fondecyt Regular #1085278, the FONDAP Center for Astrophysics #15010003 and the 
Basal CATA PFB-06. 
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Table 1. Iron abundances and AV^™^ for the 15 GCs in our sample. 



ID Alias yfe^mp AV^g^ [Fc/H] 

RHS97(ZW84)'^ RHS97(CG97)^ COQ*' KI04'= 
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^Iron abundances by Rutledge et al. (1997, RHS97) in the Zinn & West (1984, ZW84) and in the Carretta & Gratton 
(1997, CG97) metalhcity scale. 

''Iron abundances by Carretta et al (2009, C09). 

=Iron abundances by Kraft & Ivans (2003, 2004, KI04). 
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Table 2. Evolutionary predictions for stellar structures constructed at fixed cluster age 
(12 Gyr), global metallicity and assuming different helium abundances. 
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^Stellar mass of the progenitor at the tip of 
the RGB. 

'^Surface extra helium after the first dredge- 
up. 



